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Abstract 

Future IR missions should give some priority to high resolution spectroscopic 
observations of the set of far-IR transitions of OH. There are 15 far-IR lines 
arising between the lowest eight rotational levels of OH, and ISO detected 
nine of them. Furthermore, ISO found the OH lines, sometimes in emission 
and sometimes in absorption, in a wide variety of galactic and extragalactic 
objects ranging from AGB stars to molecular clouds to active galactic nuclei 
and ultra-luminous IR galaxies. The ISO/LWS Fabry-Perot resolved the 
1 19pm doublet line in a few of the strong sources. This set of OH lines 
provides a uniquely important diagnostic for many reasons: the lines span a 
wide wavelength range (28.9pm to 163.2pm); the transitions have fast 
radiative rates; the abundance of the species is relatively high; the IR 
continuum plays an important role as a pump; the contribution from shocks 
is relatively minor; and, not least, the powerful centimeter-wave radiation 
from OH allows comparison with radio and YLBI datasets. The problem is 
that the large number of sensitive free parameters, and the large optical 
depths of the strongest lines, make modeling the full set a difficult job. The 
SWAS montecarlo radiative transfer code has been used to analyze the 
ISO/LWS spectra of a number of objects with good success, including in both 
the lines and the FIR continuum; the DUSTY radiative transfer code was 
used to insure a self-consistent continuum. Other far IR lines including those 
from H 2 0, CO, and [01] are also in the code. The OH lines all show features 
which future FIR spectrometers should be able to resolve, and which will 
enable further refinements in the details of each cloud’s structure. Some 
examples are given, including the case of SI 40, for which independent 
SWAS data found evidence for bulk flows. 

1 The Far Infrared Transitions of OH 

Storey, Watson and Townes (1981) made the first far infrared detection of OH in the 
interstellar medium: the two 1 19pm lambda-doubled lines between the ground and first 
excited states, which they discovered in absorption from Sag B2, and in emission from 
the shock in Orion KL. Altogether there are 15 far infrared lines between the lowest 8 
rotational levels of OH. The lines (each involving six transitions between the hyperfine- 
split rotational levels) are at wavelengths of (approximately) 28.9pm, 34.6pm, 43.9pm, 
48.8pm„ 53.3pm, 65.1pm, 71.2pm, 79.2pm, 84.4pm, 96.3pm, 98.7pm, 115.4pm, 
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1 19.9pm, 1 34.8pm, and 163.1|um. The upper state excitation temperatures for these 
lines range from 120K to 618K. The dipole moment for OH is large, 1.668 Debye (for 
comparison the CO dipole moment is 0.1 12 Debye), and the radiative rates for OH 
transitions are generally fast. For example, the 1 19pm fundamental transition rate is 
about 0.1 sec 1 . But the FIR OH transitions also include some cross-ladder lines whose 
radiative rates are one hundred times weaker, providing a dataset of neighboring, far IR 
lines which frequently include both optically thin and very optically thick features. The 
OH analyses have an additional resource from which to draw: the strong hyperfme radio 
wavelength transitions that OH has in its ground-state, and which have been extensively 
observed. Maser and/or mega-maser activity is seen in many of the stars and galaxies 
which ISO observed. 

The general properties of interstellar OH are known from thermal OH emission studies 
done at radio wavelengths, as well as from the far IR observations. Typically OH has the 
following range of properties: N OH /N H2 = 0.1-3 xlO' 7 ; N OH L = 2-100 x 10 15 cm' 2 ; T x = 
100 - 275K; and N H2 =0.1 - 3 xlO 7 cm' 3 (e.g., Watson et al. 1985; Jones et al. 1994). In 
maser galaxies, such as the ones we observed with ISO/LWS, the OH masers regions 
typically have somewhat different properties: N oh / N H2 = 10' 7 - 10' 8 ; N OH L ~10 15 cm' 2 ; 
T x = 40K - 50K; N H2 ~ O.l-lxlO 4 cm' 3 (e.g., Henkel and Wilson, 1990). The presence of 
OH megamasers allow for VLBI observations in AGN, and results indicate that the sizes 
of the emitting regions are compact: ~ from a few to tens of parsecs. For example, in 
Arp220, one of the galaxies I discuss below, they seem to surround the AGN, with H 2 
densities ~ 10 6 cm' 3 (Lonsdale et al. 1994; Sk inn er et al. 1997). 

Four of the relatively strong far IR OH transitions involve the ground state. In warm 
clouds the molecules absorb the strong IR dust continuum, populating effectively some of 
the higher lying levels. As was conclusively shown by Sylvester et al. (1997), the 18 
centimeter radio maser emission in evolved stars is pumped by absorption of the 34pm 
dust continuum, while in AGN Skinner et a l (1997) proved the effectiveness of the IR 
pumping of mega-masers in Arp220. The presence of a strong far IR continuum affects 
all of the far IR lines, to varying degrees, and means (besides making the models more 
complicated) that the set of far IR OH lines also provide a sensitive measure of the local 
continuum conditions. One final point is worth noting about OH: unlike other some 
species commonly used as far IR probes of the interstellar medium, like CO, H 2 0, or OI, 
OH emission from PDRs or shocks does not contribute a relatively dominant amount of 
luminosity to these other processes. In the cases I model below, just the warm gas from a 
dusty molecular clouds is adequate to explain the strengths we observe. OH, therefore, 
has a powerful combination of features that makes it a very useful species for dis- 
entangling cloud properties across a very wide range of conditions. 
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2 ISO/LWS Observations of Extragalactic OH 

The ISO/LWS Extragalactic Science team has seen OH in fifteen galaxies, and has 
obtained potentially useful limits on approximately one hundred other galaxies observed 
by ISO. Our observations include a set on the 34.6pm “IR pump” transition from the OH 
2 PI 3/2 ground state, which was obtained using the ISO/SWS spectrometer. This line, and 
the less important 28.9pm line between excited upper states, are the only OH far IR lines 
not in the LWS wavelength coverage. ISO of course has also seen OH in numerous 
galactic locations, including Sag A, and Sag B2, and in particular in the evolved star 
IRC+ 10420 — the first source to demonstrate that the 34pm continuum can effectively 
pump the OH maser lines in these stars. The LWS extragalactic detections of OH are in 
the following sources: Arp 220, CenA, IRAS 17208-00 14, IRAS20 100-4 156, M82, 
Mkn 23 1 , Mkn 273, NGC 253, NGC 891, NGC 1068, NGC 1614, NGC 3690A, 

NGC 4945, NGC 7469, and 3Zw35. Figure 1 shows all the eight OH lines detected in 
Arp220, along with a sample of the lines detected in other objects to illustrates some of 
the morphological variety. 

2.1 General Characteristics of the Observed Extragalactic OH Lines 
Perhaps the most striking observation about the set of extragalactic lines measured is the 
wide morphological range of behavior they display, even though all arise in infrared 
bright galaxies with either active star formation, an active nucleus, or perhaps both. In 
Arp 220, for example, every OH line is seen in absorption except the longest wavelength, 
163pm, line which is seen weakly in emission (Fischer et al. 1998) . By contrast, NGC 
1068, another AGN, has every detected line seen in emission , even the strong 1 19pm 
transition between ground and the first excited state (Spinoglio et al. 1999). NGC 253, a 
nearby starburst, has some OH lines in emission and some in absorption (Bradford et al. 
1999), while in M82, the infrared bright, prototype starburst galaxy, the lines’ equivalent 
widths are so small continuum that even with our high signal-to-noise ratio only the 
1 19pm line has been conclusively seen at all, in absorption (Colbert et al. 1999). From 
the analysis of this diverse set of line strengths, several useful preliminary generalizations 
may be drawn for the different categories, as summarized below. 

ULIGs: In the case of Arp 220, the VLBI megamaser studies together with the strength 
of the 34pm pump absorption provide some strong physical constraints. The OH lies in 
numerous small clouds which surround the AGN, and which are pumped by the local, 
warm, far IR continuum. Through further modeling I hope to arrive at a better sense of 
the shape of that continuum: is it starburst-like, or more AGN like, and is it characteristic 
of all ULIGs? 

AGN: In the case of the bright Seyfert galaxy, NGC 1068, analyses of the strong atomic 
lines (Spinoglio et al. 1999) show the substantial presence of a PDR line emission 
component, along with the high ionization lines from the AGN component in the ISO 
beam, and a starburst component, which is seen as well in other AGN. However the 
[CIIJ/FIR ratios are strange — often less than in PDRs. The OH lines, all seen in 
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Figure 1: The 8 ISO/LWS OH lines seen in Arp 220 (the 34pm SWS-band line is not 
shown here). Under each OH line in Arp220, an OH line selected from another other 
ISO galaxy (or IRC+ 10420) is shown for comparison at the corresponding wavelength to 
illustrate some of the variety in the intensities observed. 
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emission, might help sort out the density and geometry of the clouds (for example, to see 
if the PDR regions have smaller than average filling factors) - and provide clues to their 
relationship to the active nucleus. 

Infrared Bright Galaxies: We find that the OH 119 pm ground fundamental transition is 
always in absorption in these galaxies, as is the 53pm line, while the 163pm feature is 
always seen in emission. Otherwise there does not appear to be any consistent behavior 
in the lines from different objects: in some sources they are seen in emission, in others 
they are in absorption. Saraceno et al. (1996) and Benedetti et al. (2000), among others, 
have noted there seems to be a dearth of H 2 0 emission in some galactic clouds. Also, 
the [CII] fluxes are very low, perhaps due to low gas heating (but this is not conclusively 
demonstrated). They note that the neutral oxygen [01] 63pm line is often abnormally 
weak and may be self absorbed. Finally, the SWAS satellite found that 0 2 is very weak 
or absent (Goldsmith et al. 2000). OH plays a key role in the chemistry of the ISM, is 
sensitive to the temperatures and radiation fields, and its abundance and distribution 
should help in the analyses of all these issues. 


3 Modeling of the Far Infrared OH Lines 

3.1 The Montecarlo Code for Lines: the “DUSTY” Code for Continuum 
I used the one-dimensional montecarlo radiative transfer code developed by the SWAS 
mission (Ashby et al. private comm) to model the OH line strengths. This code, a 
modification of the original Bernes code, adds a treatment of continuum photons from 
dust mixed in with the gas — a particularly essential feature for OH, which is pumped in 
many cases by absorption of 34pm continuum. In addition, the code corrects for some 
previous errors encountered at large optical depths, also an issue of importance for OH 
which has a very strong matrix element. Finally, the code includes an ability to handle a 
wide range of molecules besides rigid rotors. The montecarlo code takes as input a series 
of concentric shells, each of which is specified as to size, gas and dust temperature, H 2 
density, velocity and velocity width, and molecular abundance relative to H 2 . The model 
as a whole also assumes a (specifiable) dust emissivity. In all the modeling, the dust is 
assumed to be 1% of the gas, and to have the same temperature as the gas everywhere in 
the cloud. The montecarlo code calculates the populations of the molecular levels in 
each shell. The output is then fed into a radiative transfer code that calculates the line 
profiles as seen by an external observer looking at the cloud. 

The montecarlo results confirm the fact, noted above, that the OH lines are often optically 
thick. The peak optical depth at line center for the 1 19pm line exceeds 100 in some 
cases. As a result, many of the lines in many situations show self-absorption. As the 
model is tuned to fit the data by increasing the column density, these features can turn an 
emission line into an absorption line, one while strengthening the emission from a weaker 
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neighbor. The full set of nine line intensities, which span a range of wavelengths and 
optical depths, enable us to derive a rather detailed self-consistent picture of the cloud 
conditions solely from their intensity ratios, without the need for velocity information. 
But, as seen in the example below, there is also a wealth of information in the line shapes. 

ISO/LWS also observed lines of CO, H 2 0, and [01] in these galaxies, with varying 
degrees of success. The SWAS montecarlo code can also predict the line emission from 
these species. In general the model predictions for these lines are in overall agreement 
with the observations; the uncertainties are due to uncertainties in the ISO line fluxes 
themselves, in the assumed molecular abundances, and in the amounts of possible line 
“contamination” from shocks and PDRs in the beam. Overall the results add further 
confidence to the models. The [01] 63 pm line is clearly self absorbed in several 
instances, confirming the suspicion that the line intensity is sometimes very weak due to 
self-absorption (Saraceno et al. ) 

The montecarlo output is not strictly internally self-consistent; the input parameters need 
not conserve luminosity between shells, for example. To obtain this self consistency, and 
to guarantee that the final cloud model was also consistent with the observed far infrared 
continuum flux, I used the DUSTY code (Ivezic and Elitzur, 1997) to model the 
continuum and generate a set of shell parameters that provided this consistency. Then I 
iterated the DUSTY model with the montecarlo line outputs. While this technique does 
not give a unique solution for the cloud structure, it does provide a canonical model 
consistent with the observations. 


3.2 Some Modeling Results 

3.2. 1 S 140 - A Molecular Cloud with Bulk Inflow Motions 

Ashby et al. (2000) used the SWAS satellite to observe S140, and the SWAS montecarlo 
code to model the shape of the observed submm H 2 0 line. They are able to obtain a 
detailed, if not entirely unique, model of the cloud. From their set of models they 
conclude, for example, that the cloud radius is 0.44pc, has an inner temperature of ~70K 
an inner hydrogen density of 1.4x106 cm' 3 , a radial profile of temperature which varies 
like r'° 5 , and a density profile varying like r'° 8 . They also conclude that “significant bulk 
flow” is required to explain the H 2 0 line shape, but because of the small optical depth of 
the 557 GHz line they could not differentiate between infall and outflow. Although the 
ISO observations of SI 40 were only able to set weak limits on the OH lines (Aannestad 
and Emery, 1998 ), the detailed nature of the SWAS models made it a useful check. 
Taking the SWAS cloud parameters for the case of bulk infall motions, I used the 
montecarlo code to calculate the strengths and line shapes for the full set of OH lines; I 
calculated the far IR lines of H 2 0 and the [OI] lines as well, for comparison with the ISO 
values. The reasonable results obtained from the model of S 140, as confirmed with the 
SWAS (and ISO) observations, provide some confidence in the galaxy modeling. Figure 
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2 shows the results for three of the OH lines, on the same scale: the 163 jam line (top), the 
79pm line (middle), and the 1 19pm line (bottom). It is clear that high spectral resolution 
observations of the far IR water lines can readily distinguish infall from outflow, because 
these lines are optically thicker than the submillimeter line. 

3.2.2 Arp220 - The Molecular Cloud Component of a Peculiar Ultra Luminous Galaxy 
Arp220 is a particular challenge, because so many OH lines are seen, and every one of 
them (except the line at 163 pm) is seen in absorption — something that happens in no 
other known object. Arp220 is also unusual in general in that its spectrum is charac- 
terized by numerous molecular absorptions; even the [01] 63pm fine structure line is seen 
in absorption, while the [CII] 157pm line is found only in weak emission (e.g., Fischer et 
al. 1998). Skinner et al. (1997) showed that the 34pm OH absorption could pump the 
OH megamasers in Arp220, and their analysis supported the model in which many small, 
dense molecular clouds circulate around the nucleus. Suter et al. (1998) attempted to 
model the OH line absorptions with some less complex radiative transfer models, and 
were driven to consider some unusual non-LTE situations to explain the observations; . 

I used the SWAS montecarlo code successfully to model (to first order) all of the 
observed OH lines in Arp220. Figure 3 shows a sample of these model line shapes. In 
order to get absorption in all the lines, including the 79pm line to the ground state 2 PI 3/2 
from the lowest level of the 2 PI 1/2 ladder, while still getting weak emisson at 163 pm, it 
was necessary to have falling within the ISO/LWS beam a combination small dense 
clouds, and a few large, giant molecular clouds. Iterating with the DUSTY code provides 
a way to obtain a continuum shape that fits the observed infrared continuum, though there 
may also be a continuum component in our beam arising from clouds without much OH. 
It was not necessary in this process to resort to the non-LTE scenario postulated in Suter 
et al. 

3.3.3 NGC 253 - The OH Emission from a Nearby Infrared Bright Galaxy 

In NGC 253, the two strongest OH lines are seen as emission lines: the 79pm and 163pm 
features, each of which is about twice as strong as the 1 19pm fundamental absorption 
(Bradford et al. 1999). Only two other OH lines are also seen in this galaxy, and as a 
result the model’s constraints are weaker than for Arp220. The montecarlo modeling of 
these OH lines implies that a few giant molecular clouds, with r ~ lOOpc, can explain the 
observations. The clouds have H 2 densities averaging approximately 3xl0 3 cm' 3 and an 
OH abundance relative to H of about 5xl0' 8 . ISO/LWS also observed the 1 19pm line in 
NGC 253 with the high resolution Fabry-Perot, which Bradford (2001) has successfully 
modeled. 

3.3.4 NGC 1068 - OH Emission from AGN 

NGC 1068 is a remarkable extragalactic source in its OH spectrum, because it is the only 
galaxy observed in which all the detected transitions are found in emission. The code 
accurately predicts the observed line flux ratios to 40%. Based on the montecarlo 
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Figure 2: The hydrogen density distribution in a nominal giant molecular cloud in the 
galaxy NGC 253. This distribution was used in the montecarlo code modeling. The 
inner 50 parsecs of the GMC has warm, 250K, gas and dust, at constant density, while the 
bulk of the cloud has much cooler, 35K material. 



Figure 3: Predicted OH line shapes from infalling gas in SI 40, based on cloud profiles 
derived from SWAS data. The figure shows the 163 pm line (top), the 79 pm line 
(middle), and the 119 pm line (bottom). While all three lines show evidence for infall 
from their highly asymmetric shapes, only the latter two show characteristic absorption: 
in the redshifted (near side) material. The velocity range across the figure is ±10 km/sec, 
reflecting the fact that the maximum infall velocity in the model (and seen by SWAS) 
was 7.2 km/sec. 
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Figure 4: Montecarlo modeling results of three OH lines in Arp 220. The OH lines at 
1 19pm, 79 pm, and 53 pm are displayed. All the OH lines are seen in absorption in Arp 
220, except the weakly emitting 163 pm line - the only known object with this property. 
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modeling, it seems the clouds that produce these lines in NGC1068 must be relatively 
small and dense, and heated from the inside. They are about 0.2pc in radius, with 
densities of ~ 10 4 cm' 3 at the outer edges, increasing towards the center with a power law 
behavior of (R/Ro) 1 25 . The temperature in the outer shells is about 25K, increasing 
inwards approximately with a power law dependence (R/Ro) 1 47 . The total column 
density of H 2 in each cloud is 1.5xl0 24 cm' 2 , and the relative abundance of OH in the 
model is approximately 10' 7 . OH is taken to be absent in the hotter portions of the cloud, 
where T>300K. With these high column densities the strong lines of OH are very 
optically thick, and precise radiative transfer calculations like the present ones are 
absolutely essential. At a distance of 16.2Mpc, the total number of equivalent number of 
clouds needed to produce the observed absolute flux in the lines and the continuum is 
~3xl0 7 . The limits to the observed fluxes in all the weaker lines are consistent with this 
modeling. It is worth noting in this context that in some galaxies the weaker OH lines in 
NGC 1068, like the 53pm line, are amongst the strongest, but the montecarlo modeling 
can successfully account for these differences. The montecarlo code also successfully 
reproduces the continuum emission seen from NGC 1068 to within a factor of 2-3 in 
absolute flux density across the entire LWS spectrum. The total mass in such an ensemble 
of clouds in NGC 1068 is approximately 9xl0 9 M®. 


4 Conclusions to Date 

What general conclusions might we hope to draw from this large set of ISO OH 
observations? Although the analysis is still underway, it appears that the infrared 
luminous galaxies can be grouped into three general categories based on the relative 
strengths of their infrared OH lines: (1) relatively normal galaxies like NGC 253 and 
M82, from which the OH lines are seen in both emission and absorption, and which are 
dominated by starburst activity in giant molecular clouds; (2) NGC 1068 and other AGN, 
whose OH is seen in emission, and whose clouds must be quite small and dense 
(hydrogen densities up to 10 s cm' 3 in the cloud cores); and (3) Arp220, and perhaps other 
ULIGs with active starbursts as well as an active nucleus, for which a combination of a 
few giant molecular clouds and a collection of moderately dense cloudlets are required. 
(The ISO/LWS spectrum of Mkn 231 is much noisier than Arp 220's, but hints at similar 
behavior in its OH.) While I have made substantial progress using only the intensities of 
these lines, velocity resolved spectra will enable even further refinements. Future far 
infrared missions with sensitive, spectroscopic capabilities should find that the set of OH 
lines can provide wealth of information needed to unravel the structures of complex 
clouds. 


Acknowledgments : 

The modeling work reported here was done in close collaboration with the SWAS 
montecarlo team, especially Matt Ashby and Gary Melnick. The ISO observations and 
their analysis represents a large effort by the ISO/LWS extragalactic science team, in 


107 


Second Workshop on New Concepts for Far-Infrared and Submillimeter Space Astronomy 


particular Jackie Fischer, Chris Dudley, Shobita Satyapal, Luigi Spinoglio, Gordon 
Stacey, Matt Bradford, Matt Malkan, and Sarah Higdon. I particularly want to thank the 
Istituto di Fisica dello Spazio Interplanetario (IFSI) in Rome, Italy, and Paolo Saraceno of 
the host Infrared Group, for their gracious hospitality during a summer visit during which 
time I did much of this analysis. The research was sponsored in part by NASA Grant 
NAG5- 10659. 

References: 

Aannestad, P.A. and Emery, R, 1998, BAAS, 1 92, 1009. 

Ashby, M., et a!., 2000, Ap.J. (Letters), 539, LI 19. 

Benedettini, M., Giannini, T., Nisini, B., Tommasi, E., Lorenzetti, D., Di Giorgio, A. M., 
Saraceno, P., Smith, H. A., and White, G. J., 2000, A&A , 359^ 148. 

Bradford, C. M., Stacey, G. J., Fischer, J., Smith, H. A., Cohen, R. J., Greenhouse, M. 

A., Lord, S. D., Lutz, D., Maiolino, R., Malkan, M. A., and Rieu, N. Q., 1999, 

The Universe as Seen by ISO. Eds. P. Cox & M. F. Kessler. ESA-SP 427. 

Bradford, C.M., 2001, Ph.D.thesis. 

Colbert, James W., Malkan, Matthew A., Clegg, Peter E., Cox, Pierre, Fischer, 

Jacqueline, Lord, Steven D., Luhman, Michael, Satyapal, Shobita, Smith, Howard 
A., 1999, Ap,J., 51 1, 721. 

Fischer, J., Satyapal, S., Luhman, M., Melnick, G., Cox, P., Cemicharo, P., Stacey, G., 
Smith, H.A., Lord, S., and Greenhouse, M., 1998, in Proceedings of the “ISO to 
the Peaks Workshop: The First ISO workshop on analytical Spectroscopy”, 
Madrid, Spain. 

Goldsmith, P., et a/., 2000, Ap.J. (Letters), 539, LI 23. 

Henkel, C., and Wilson, T.L., 1990, A&A, 22! 91 431. 

Ivezic, Z., and Elitzur,M., 1997, M.N.R.A.S., 287, 799. 

Jones, K.N., Field, D., Gray, M.D., and Walker, R.N.F., 1994, A&A, 288 581. 

Lonsdale, C.J., Diamond, P.J., Smith, H.E., Lonsdale, C.J., 1994, Nature, 37Q 117. 
Saraceno, P., et al. 1996, A&A, 315 L293-L296. 

Skinner, C.J., Smith, H.A., Sturm, E., Barlow, M.J., Cohen, R.J., and Stacey, G.J., 1997, 
Nature, 388 472. 

Spinoglio, Luigi, Stacey, Gordon, and Unger, Sarah J, 1999, Ap.J., 51 1, 721. 

Spinoglio, L., Suter, J. P., Malkan, M. A., Clegg, P. E., Fischer, J., Greenhouse, M. A., 
Smith, H. A., Pequignot, D., Stacey, G. J., Rieu, N. Q., Unger, S., J., and 
Strelnitski, V., 1999, The Universe as Seen by ISO. Eds. P. Cox & M. F. Kessler. 
ESA-SP 427. 

Storey, J.W.V., Watson, D.M., and Townes, C.H., 1981 , Ap.J. (Letters), 244, L27. 

Suter, J. P., Harvey, V., Cox, P., Fischer, J., Lord, S. D, Melnick, G. J, Satyapal, S., 

Smith, H. A., Stacey, G., Strelnitski, V., and Unger, S. J., 1998, BAAS, 193, 9002. 
Watson, D.M., Genzel, R., Townes, C.H., and Storey, J.W.V., 1985, Ap.J, 298 316. 


108 


